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ABSTRACT 

I use high-quality HST spectrophotometry to analyze the calibration of three 
popular optical photometry systems: Tycho-2 S^Vt, Stromgren uvhy, and John- 
son UBV. For Tycho-2, I revisit the analysis of Mafz Apellaniz (2005a) to 
include the new recalibration of grating/aperture corrections, vignetting, and 
charge transfer inefficiency effects produced by the STIS group and to consider 
the consequences of both random and systematic uncertainties. The new results 
reaffirm the good quality of both the Tycho-2 photometry and the HST spec- 
trophotometry, but yield a slightly different value for ZP^t-Vt of 0.033 ± 0.001 
(random) ±0.005 (systematic) magnitudes. For the Stromgren f , 6, and y filters 
I find that the published sensitivity curves are consistent with the available pho- 
tometry and spectrophotometry and I derive new values for the associated ZP^.j^ 
and ZPmi . The same conclusion is drawn for the Johnson B and V filters and the 
associated ZP^.y. The situation is different for the Stromgren u and the Johnson 
U filters. There I find that the published sensitivity curves yield results that are 
inconsistent with the available photometry and spectrophotometry, likely caused 
by an incorrect treatment of atmospheric effects in the short-wavelength end. I 
reanalyze the data to produce new average sensitivity curves for those two filters 
and new values for ZPc^ and ZP u_b- The new computation of synthetic U — B 
and B — V colors uses a single B sensitivity curve, which eliminates the previous 
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unphysical existence of different definitions for each color. Finally, I find that if 
one uses values from the literature where uncertainties are not given, reasonable 
estimates for these are 1 — 2% for Stromgren b — y, mi, and Ci and 2 — 3% for 
Johnson B — V and U — B. The use of the results in this article should lead to a 
significant reduction of systematic errors when comparing synthetic photometry 
models with real colors and indices. 

Subject headings: space vehicles: instruments — stars: fundamental parameters 
— techniques: photometric — techniques: spectroscopic 

1. Introduction 

The data explosion of the last decade has produced large quantities of photometric 
measurements, either from dedicated ground-based projects (e.g. 2MASS, Skrutskie et al. 
1997), compilations from different sources (e.g. GCPD, Mermilliod et al. 1997), or dedicated 
space-based missions (e.g. Hipparcos, ESA 1997). The future promises an expansion of the 
phenomenon with the final publication of data from projects such as the Sloan Digital Sky 
Survey (York et al. 2000) or GALEX (Bianchi et al. 1999) and new missions such as GAIA 
(Mignard 2005). 

The photometry from space-based missions is obtained with a fixed instrumental setup, 
does not suffer from perturbing atmospheric effects, and can be supported by a dedicated 
calibration program. On the other hand, observing from space also has its negative effects 
(e.g. the existence of Charge Transfer Inefficiency, or CTl, in CCD detectors due to radi- 
ation damage, Kimble et al. 2000) but, overall, the available data from recent space-based 
missions have very good photometric characteristics. The dedicated ground-based projects 
of the last decade also benefit from the uniformity of the instrumental setup, observing site, 
and reduction techniques and, though affected by the atmosphere, they can implement ex- 
tensive calibration programs and quality control mechanisms to ensure the stability of the 
obtained data (see, e.g. Cohen et al. 2003, Smith et al. 2002). Large dedicated photometric 
projects (either space- or ground-based) also benefit from the volume of the datasets and 
their coverage of a significant part of the sky (or all of it), which allow for statistical tests 
to detect possible systematic effects in the data. 

The situation is quite different for compilations of older ground-based data (see, e.g. 
Lanz 1986; Hauck & Mermilliod 1998), which typically include tens or hundreds of studies, 
each one of them with tens or hundreds of objects. The large number of sources implies 
different detectors, filters, observatories, and reduction techniques, which undoubtedly in- 
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troduce a scatter in the data and, quite possibly, systematic errors (see Stetson 2005 and 
references therein for a detailed discussion on the problems associated with the reduction 
of ground-based photometry). Furthermore, in many cases the photometry obtained in the 
"standard" systems (e.g. Johnson or Stromgren) was established using equipment with lower 
precision than the one available today and using standard stars with a limited range of col- 
ors. Therefore, it is possible that some systematic errors may have been introduced from 
the very beginning, thus comphcating a recalibration of the data. Those effects led Bessell 
et al. (1998) to conclude that (a) standard systems may not represent any real linear sys- 
tem and (b) we should not be reluctant to include ad-hoc corrections of a few percent to 
achieve an agreement between the observed and the synthetic photometry in such systems. 
This rather pessimistic view may lead somebody to ask oneself why should we bother about 
the recalibration of photometric systems half-a-century old. The answer to that question is 
that many of the current calibrations are eventually tied up to those standard systems, as 
evidenced by the fact that the most common way to simply characterize the photometric 
properties of an object in a research proposal is to give its Johnson V magnitude and B — V 
color. If we renounce to calibrate standard systems accurately we are condemned to have 
systematic errors propagating down our reduction procedure. 

In this paper I take a more optimistic view of the problem by stating two points. In the 
first place, even though it is true that some standard systems cannot be strictly characterized 
by a sensitivity (or throughput) curve and a zero point for each magnitude and/or color, 
it is also true that does not stop us from trying to derive the corresponding function and 
value that minimize the scatter of the data. In the second place, some of the calibration 
problems encountered in the past may have been due not to the photometry itself but to the 
data used to calibrate it, e.g. systematic errors in the measured spectrophotometry or in the 
assumed stellar model parameters. Therefore, the use of more modern calibration data may 
increase the accuracy. In other words, I believe it is possible to eliminate at least some of the 
systematic errors in the photometric calibration and to significantly lower the uncertainties 
in the zero points. 

In paper I (Maiz Apellaniz 2005a), I used HST spectrophotometry to analyze Tycho-2 
photometry in order to check the accuracy of its published sensitivity curves and to calculate 
the zero points. In this paper, I start by revisiting those data and then I use similar techniques 
to analyze the two most common optical standard systems: Stromgren uvby and Johnson 
UBV. In all cases, Vega will be used as the reference spectrum to calculate magnitudes, 
colors, and indices (see Appendices). 
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2. Spectrophotometry 

A straightforward method to test the sensitivity curves and to calculate the zero points 
of filter systems is to observe a number of stars with well-characterized and sufficiently 
different spectral energy distributions (SEDs) and to compare the measured photometric 
magnitudes/colors/indices with those derived from synthetic photometry^. If the sensitivity 
curves are correct, a plot of the difference between the measured and synthetic values as a 
function of color or magnitude should yield a straight line with zero slope and an intercept 
equal to the zero point (see Appendices). If a slope is present in such a plot, then the assumed 
sensitivity curves may be incorrect due to e.g. an inaccurate characterization of atmospheric 
effects (which enter in the sensitivity curve for ground-based observations in the reduction 
process to extrapolate to zero air masses) or an imprecise calibration of the detector. Of 
course, it is also possible that the cause of a non-zero slope lies in the SED library, and not 
in the photometry itself. Therefore, a necessary preliminary step in the process is to ensure 
that the SEDs are as accurate as possible. 

The Space Telescope Imaging Spectrograph (STIS) onboard Hubble Space Telescope 
(HST) was able to provide during its seven years of operation the most accurate to-date 
spectrophotometry in the 1 150-10 200 A range. The accuracy of its absolute flux calibration 
in the V band is 1.1%, of which 0.7% corresponds to the value of the absolute flux of Vega 
at the V band (Megessier 1995) and 0.8% to uncertainties^ in the photometry (Bohlin 2000; 
Bohlin & Gilliland 2004b). On the other hand, the accuracy of its relative flux calibration 
with respect to the V band depends on the uncertainty of the temperature of the three 
white dwarfs used as primary calibrators. This can be as high as 1-3% in the FUV but 
in the optical it is significantly lower due to the degeneracy of the spectral slopes at high 
temperatures. The accuracy of an individual STIS spectrophotometric observation is also 
hmited by the instrument repeatabihty, which is 0.2-0.4% for the wavelength range in which 
I am interested in this paper (Bohlin 2000; Bohlin & Gilliland 2004a). 

In this article I will use the same two STIS samples as in Paper I: the Next Generation 
Spectral Library (NGSL; Gregg et al. 2004) and the Bohlin sample (Bohlin et al. 2001; Bohlin 
& GiUiland 2004a). The reader is referred to Paper I for details. The NGSL sample is large 



^Alternatively, one can substitute the observed SEDs for synthetic (or model) ones, but doing so raises 
two issues: (a) the possible existence of errors in the SED modehng and/or in the extinction characterization 
and (b) the need to obtain additional data to measure the relative flux calibration between the model SEDs 
and the real reference SED in Eq. Al in order to derive not only color zero points but also the zero point for 
a reference magnitude such as V. For those reasons, observed SEDs are in general preferred to model ones. 

^All of the uncertainties in this article are quoted as 1 a. 
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(> 250 non-variable stars), is made up of relatively bright stars {V between 2.5 and 11.5), 
and includes a large variety of SED types (sampling diverse temperatures, metallicities, and 
gravities), but was obtained with a STIS configuration (52X0.2 long slit at the El detector 
position and reading only a subarray of the CCD) which introduces some complications in 
the calibration. The Bohlin sample is smaller (19 stars), consists of dimmer objects (only 
2 are brighter than V = 9, implying that only some of the stars have reliable multicolor 
photometry available in the standard systems), and is made up mostly of early- type stars 
(only 5 stars have B — V> 0.0). The latter sample, however, uses the preferred STIS 
setup for absolute photometry (52X2 long slit at the center of the detector and reading 
the full CCD) and, in most cases, have repeated observations, thus allowing for a more 
straightforward calibration. 

Some of the calibration issues regarding the NGSL sample were discussed and dealt with 
in Paper I. Since that article was pubhshed, the STIS group (Proffitt 2005) has identified a 
number of issues that prompted me to reprocess the data in order to improve the calibration. 
On the first place, it was discovered that the effective throughput as a function of wavelength 
for a given long slit (e.g. 52X2 or 52X0.2) was not the same for all STIS gratings. This effect 
has been corrected in the STIS pipeline by introducing a grating/ aperture correction table 
and making the appropriate changes in the calstis software. Second, a new L-flat that 
incorporates the effects of vignetting and that should improve the flux cahbration at the El 
detector position has been calculated. Finally, an error in the algorithm that calculates the 
CTI for subarrays has been fixed. All of the above issues have no effect on the broadband 
colors of the Bohhn sample but can produce changes at the 1-2% level for the broadband 
colors of the NGSL sample. For that reason, before proceeding with an analysis of the 
Stromgren and Johnson systems, I study the effects of the new calibration on the paper I 
results for Tycho-2 photometry. 



3. Tycho-2 BtVt 

As previously mentioned, the purpose of paper I was to test the pubhshed sensitivity 

curves for the two Tycho-2 filters, i?T and Vr (Bcssell 2000), and to calculate their zero 
points. I found out in that paper that the sensitivity curves were indeed correct and derived 
a value of ZFb^ = 0.078 ± 0.009 from the Bohlin sample and of ZFb^-Vt = 0-020 ± 0.001 
from the NGSL sample. The uncertainties quoted there arc only the random ones and were 
derived by inverse variance weighting from the photometry. The new calibration of the NGSL 
data can (and indeed does, as we will see next) affect the value for ZPbj—Vt but not that of 
ZPbt, which was derived from the Bohlin sample. 
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I have reprocessed the NGSL data through calstis and retained only the non-variable 
stars with a 5 600 A jump (the point where the G430L and G750L gratings overlap) of 
less than 3% and with uncertainties in the measured i?T and Vr magnitudes smaller than 
0.06 mag. This left 255 objects (as opposed to 256 in Paper I) for which the synthetic (or 
spectrophotometric) and Vt magnitudes were computed. 1 show the difference between 
the photometric values (from the Tycho-2 catalog) and spectrophotometric values (from the 
STIS spectra) of — Vt in Fig. 1, which is the equivalent to Fig. 2 in Paper I^. 

The main conclusion remains the same as in the previous paper: there is a very good cor- 
respondence between the photometric and spectrophotometric B-y — Vr colors, as evidenced 
by the nearly-flat distribution of points in Fig. 1. The normalized histogram in Fig. 2 also 
corroborates that assertion: the residuals appear to be symmetrically distributed around the 
zero point and the standard deviation of the histogram is 1.04, as in Paper I. Here I have 
also performed an additional test by fitting a weighted straight line to the data in Fig. 2, 
which is shown in green. The slope of the fit is nearly flat but not exactly so: 0.005 ± 0.002. 
This implies that cither the Tycho-2 sensitivities or the calibration of the NGSL data are 
not perfect. However, the effect is very small, since the deviations at the edges of the used 
color spectrum (which correspond to low-reddening O and M stars, respectively) amount to 
only ±0.006 magnitudes, which is 1/3 of the best individual photometric uncertainty for the 
colors in our sample. Also, 0.002/0.005 = 2.5, i.e. the slope is only 2.5 sigmas away from 
zero. Therefore, I do not consider necessary to derive new sensitivity curves for the Tycho- 
2 filters. Instead, this effect will be included in the error analysis below as a systematic 
uncertainty^. 

The new value I obtain for the St — Vt zero point using inverse variance weighting is 
ZP Bt-Vt — 0.033 ± 0.001, which is significantly larger than the value derived in Paper I 
(0.020 ± 0.001). The 1.3% difference between the two is within the expected range of the 
changes introduced in the new STIS cahbration (Proffitt 2005). It is interesting to note that if 
I only use the four Bohlin stars in Fig. 1 to calculate ZP Bt-Vt, I obtain a value of 0.046±0.013, 



"^The figures in this paper that plot a photometric quantity against the difference between a photometric 
and a spectrophotometric value assume an associated zero point of 0.0 in Eq. Al because their purpose is 
to calculate the real zero point by analyzing the plots themselves. 

*I consider an uncertainty in the zero point to be systematic if it arises from errors in our knowledge of 
the characteristics of our photometric or spectrophotometric setup or reduction procedure, e.g. an incorrect 
mean sensitivity curve or offsets in the flux calibration of the spectrophotometry. I consider an uncertainty 
to be random if it would still bo present in the absence of any systematic effects, e.g. effects caused by a 
finite S/N of the photometry or spectrophotometry or random variations from the mean of the sensitivity 
curve due to changing atmospheric conditions between different observations. 
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which is at a distance of 1 sigma from the value proposed here but 2 sigmas away from the 
previous one. Those values are consistent with an improvement in the new STIS calibration. 
The uncertainty of 0.001 is a purely random one but, as we have seen before, there are 
also systematic effects present. The slope to the fit described in the previous paragraph 
translates into a 1-sigma uncertainty of 0.003 magnitudes. Furthermore, the uncertainty 
of 3 000 K in the temperature of the white dwarf calibrators quoted by Bohlin & Gilliland 
(2004a) translates into an additional (systematic) uncertainty of 0.004 magnitudes, which 
was calculated using synthetic photometry of Kurucz models. Therefore, the proposed final 
value is ZP Bt-Vt — 0.033 ± 0.001 (random) ±0.005 (systematic) magnitudes. 

As for the absolute calibration of Tycho-2 photometry, 1 defer its analysis until the 
Johnson UBV photometry is analyzed. 



4. Stromgren uvby 

4.1. Description and sample selection 

The Stromgren standard system (Matsushima 1969) has been the most commonly used 
intermediate-band photometric optical system in the last 40 years. It consists of 4 filters in 
the 3000-6 000 A range, u, v, b, and y and the quantities that are presented in most works 
are: 



b-y 
nil — V — 2b + y 
ci = u — 2v + b 



(1) 
(2) 
(3) 



b — y is a color similar (but not identical) to Johnson — V or to Tycho-2 Bq^ — Vt, 
mi is an index that is most sensitive to the metallicity of a star, and Ci is an index that 
measures the strength of the Balmer jump (Stromgren 1966). Hauck & Mermilliod (1998) 
have published a catalog of photometric measurements in the Stromgren system compiled 
from photoelectric data in the literature. In its current (summer 2005) version, available 
from the General Catalogue of Photometric Data (GCPD^), it contains 66135 stars from 
572 different sources. 



^http : //obswww . unige . ch/ gcpd/ gcpd . html 
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I want to test whether the Matsushima (1969) sensitivity curves are consistent with the 
photometry in the hterature (in the sense of not having strong color terms when compared 
with the STIS spectrophotometry) and to calculate the corresponding zero points for b — y, 
mi, and Ci. Note that the Stromgren GCPD data includes only one color and two indices, 
but not magnitudes (though Johnson V is usually also given in many cases), so, as opposed 
to the Tycho-2 or Johnson cases, I only have to deal here with color/index zero points. 
We also need to consider that, as opposed to the Tycho-2 case, some of the scatter in the 
pubhshed photometric data is unavoidable due to the diversity of the sources. Therefore, 
our aims should be to obtain mean color/index values for the stars in our samples that are 
representative of that diversity and to measure the intrinsic scatter. On the other hand, 
I should not include in the sample highly discrepant values that are caused by incorrect 
reductions or misidentifications because that would misrepresent the scatter that is really 
unavoidable and may introduce biases in the interpretation of correct data. The elaboration 
of a photometric sample that is appropriate for a comparison with STIS spectrophotometry 
must allow both of those principles to be followed by selecting stars for which there are 
enough data to discriminate which values to include. With those ideas in mind, I followed 
these steps to create the Stromgren photometric sample: 

1. I started by searching in the GCPD for the original Stromgren photometry (i.e. not 
the mean values) for all non-variable stars in the NGSL sample and for all stars in the 
Bohhn sample. As I did for the Tycho-2 case, for the NGSL sample only those stars 
with a 5 600 A jump in the spectrophotometry of less than 3% were included. 

2. I calculated a weighted mean and a dispersion for b — y, mi, and ci for each star, which 
are taken to be the first estimate for the value and random uncertainty, respectively, 
for each of those color /indices. 

3. Highly discrepant individual data points (those more than 3 sigmas from the mean) 
were excluded and the mean and dispersion recalculated for b — y, mi, and Ci for each 
star. 

4. Only stars that had (a) at least three data points for each of b — y, mi, and Ci after 
the previous step and (b) random uncertainties of less than 0.035 magnitudes were 
kept in our sample and the rest were excluded. This step provides a more meaningful 
value to the statistical weight that will be derived from the random uncertainty (i.e. 
it eliminates values that are likely to be biased or of lower quality) . 

5. For the stars with Tycho-2 photometry, 1 plotted b — y vs. Bt — Vr (which should 
follow a quite tight correlation) and discarded the discrepant cases. 
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6. For each of the three indices, the remaining stars were divided into two groups: (a) 
those with 6 or more data points and (b) those with 3 to 5 data points. For the first 
group, those objects that had an uncertainty less than a critical value ex (where X = 
b — y, nil, or Ci) had that uncertainty replaced by ex- For the second group, the same 
procedure was repeated but using 2ex as the replacement in order to account for the 
lower statistical weight that should be assigned to stars with a smaller number of data 
points. There are two reasons for doing this. First, I avoid giving excessive statistical 
weights to those cases where due to roundoff or special circumstances (e.g. most 
data points coming from an identical instrumental setup and observing conditions) 
the dispersion of the data points is too low. Second, I can vary ex iteratively in 
order to get an accurate measurement of the true intrinsic scatter by demanding that, 
once a zero point has been calculated, the difference between the photometric and 
spectrophotometric values has a reduced of 1-0 (oi'i alternatively, since the number 
of stars in our sample is large, that the normalized histogram (Xphot — -^spcc)/o'x has 
a mean of zero and a standard deviation of one). The initial values of eb-y, emi, and 
£ci were obtained from Hauck & MermiUiod (1998)^. 

After the selection process was completed, I ended up with a total of 104 stars, 100 
from the NGSL sample and 4 from the Bohhn sample. In the next subsections I describe 
the results for b — y, mi, and Ci, respectively. 

4.2. b-y 

The results for b — y can be seen in Fig. 3, where I have represented the difference 
between the photometric and spectrophotometric values as a function of the photometric 
ones, in a manner analogous to that of the previous section for Bq^ — Vt. eb-y was iteratively 
determined to be 0.008 mag using the criterion stated in the previous subsection. The 
data does not manifest any significant trend with color and this shows in the slope of the 
weighted linear fit, which is —0.004 ± 0.005 mag (less than 1 sigma away from zero). Also, 
the slope of the fit has the opposite sign to the one 1 found for — Vt, suggesting that the 
(small) discrepancies lie in the photometry, not in the STIS spectra. Furthermore, no obvious 
different trend is seen between the objects in the NGSL and Bohlin samples. Therefore, I 
can conclude that the Matsushima (1969) sensitivity curves for b and y provide an accurate 
description for those filters. 



^Note that their definition of ex is slightly different from the one here. This should not be an issue 
because I only use those values as initial guesses. 
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As I did for Bt — Vr, the zero point was calculated using inverse variance weighting 
to obtain ZP^.^ = 0.007 ± 0.001. For the systematic uncertainties there are two sources, 
which I can calculate using the same techniques as for Bq^ — V^. From the weighted linear 
fit of the data there is a 1-sigma uncertainty of 0.001 mag and from the uncertainty of 3 000 
K in the white dwarf calibration I get 0.003 mag. Combining them in quadrature I obtain 
ZPb-y = 0.007 ±0.001 (random) ±0.003 (systematic) magnitudes. 

4.3. mi 

The results for mi can be seen in Fig. 3, where 1 have chosen to represent in the x axis 
the photometric b — y instead of mi. The reason for that choice is that 6 — ?/ is an almost 
monotonic function of the spectral slope in the optical range, which can be more easily used to 
detect discrepancies between the published and the real sensitivity curves, mi, on the other 
hand, is barely sensitive to spectral slope (as it should be, since the physical variable that 
produces most of its variation is metallicity) . was iteratively determined to be 0.007 
mag. The data shows a slight dependence with color and the slope of the corresponding 
weighted hnear fit is 0.010 ± 0.005 mag. Since I have already tested the validity of the 
sensitivity curves for b and y, this result indicates that the sensitivity curve for v (the other 
magnitude involved in mi) may be slightly incorrect. Note, however, that the slope is only 
2 sigmas away from zero, that the variation over the total useful b — y range amounts only 
to ±0.006 magnitudes, and that no significant differences are seen between the NGSL and 
Bohlin samples. Therefore, the situation is similar to what happens for B^ — Vr and I 
decided not to calculate any modifications in the v sensitivity curve but to include the effect 
as a systematic uncertainty. 

The zero point was calculated using inverse variance weighting to obtain ZP^^ = 
0.154 ±0.001 mag. The difference between the zero point and the linear fits yields a 1-sigma 
systematic uncertainty of 0.003 mag and the uncertainty in the white dwarf calibration con- 
tributes with less than 0.001 magnitudes, as expected from the weak dependence of mi on 
temperature. The final result is then ZP^^ = 0.154 ± 0.001 (random) ±0.003 (systematic) 
magnitudes. 

4.4. Cl 

I performed for ci an analysis similar to the one for mi using again the Matsushima 
(1969) sensitivity curves. The (fc — y)phot vs. Ci^phot — Ci^gpec plot in Fig. 5 clearly shows that a 
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color term is present, with the bluer stars above the mean and the redder ones below it. The 
measured slope is —0.045 ±0.006, which is 7.5 sigmas away from zero^ and the variation over 
the full range is ±0.027 magnitudes. Given the results I found for b — y and mi, I consider 
that variation too large to be tolerated as a systematic uncertainty. The likely culprit is one 
of the sensitivity curves and, given that 1 have previously determined that the Matsushima 
(1969) b filter curve is correct and that the corresponding v curve introduces only small 
errors in the synthetic magnitudes, I conclude that it is the u curve definition the one that 
needs to be recalculated. 

A new Stromgrcn u sensitivity curve was derived by minimization of Ci^phot — Ci^spec- 
I used a custom-made IDL code that includes the curve-fitting package written by Craig 
Markwardt^, which allows for the fitting of an arbitrary function using x^-minimization with 
restrictions on the parameters. I built the function in the following manner: (a) I selected 10 
pivot wavelengths at 75 A intervals between 3 150 A and 3 825 A. (b) The sensitivity at both 
extremes was fixed to be zero, given that we expect the result to be only slightly different 
from the Matsushima (1969) curve, (c) I placed additional constraints on the parameters to 
allow for a single-peaked function for the same reason, (d) This left nine free parameters to be 
fit, the sensitivity at the eight intermediate pivot wavelengths and ZP^i- The initial guesses 
for the parameters were derived from the Matsushima (1969) u curve, (e) The sensitivity 
curve between the pivot wavelengths was calculated using linear interpolation followed by 
smoothing with a box filter. 

An initial run of the code produced a first estimate for both the sensitivity curve and 
ZPrni, but the reduced was too large due to the underestimation of £ci- The code was 
run several times, changing until a value of 1.0 was reached for the reduced x^. The 
sensitivity curve itself was found to change little between iterations, was determined to 
be 0.009 mag and ZP^ (which here is one of the free parameters in the sensitivity-curve 
fitting procedure) to be 1.092 ± 0.002 magnitudes. The old and new sensitivity curves for 
Stromgren u are shown in Fig. 6, along with the sensitivity curves for Stromgren v and b 
and three selected spectra. The two curves for u are found to be quite similar, with the only 
exception of the short- wavelength edge, which is redder by a; 35 A for the new result. 

The new version of the (6 — |/)phot vs. ci phot — Ci^gpec plot is shown in Fig. 7. As expected, 
the new curve does a much better job of fitting the data. The color term is now absent, as a 
weighted linear fit yields a slope of —0.004 ±0.006, less than 1 sigma away from zero. Again, 
as for 6 — y and mi, no significant differences are perceived between the NGSL and Bohlin 



^The uncertainty in the slope depends on the exact value of Ea , which is determined later, 
^http : //cow . physics . wise . edu/ craigm/ idl/ idl . html 
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samples. The only outstanding issue is that the four bluest stars (two from the NGSL and 
two from the Bohlin samples) are still 1-3 sigmas above the zero point. In those cases the new 
sensitivity curve is a clear improvement (previously the same stars were 3-7 sigmas above 
the zero point) but the existence of such a residual systematic effect is likely a manifestation 
of the inherent difficulty of doing ground-based photometry to the left of the Balmer jump. 

I can conclude that the new sensitivity curve does a better job of fitting the data than 
the old one (albeit not a perfect one). As for the error budget for the zero point, the 
slope of the fit contributes 0.001 mag to the systematic uncertainty and the uncertainty 
in the temperature scale (again, calculated from the different colors of high-gravity Kurucz 
models) adds 0.004 mag. Hence, I obtain ZPcj = 1.092±0.002 (random) ±0.004 (systematic) 
magnitudes. 

5. Johnson UBV 

5.1. Description and sample selection 

The Johnson UBV standard system (Johnson & Morgan 1953) is the most commonly 
used broad-band optical photometry system. It covers a range of wavelengths similar to the 
Stromgren system but with only three filters. The quantities that are presented in most 
references are the V magnitude and the two B — V and U — B colors. Both colors are a 
function mostly of temperature and extinction but metallicity and gravity effects are also 
present, especially at the lower stellar temperature end. B — V rather monotonic function 
of temperature but for Teg > 20 000 K is almost degenerate. U — B is more appropriate to 
measure temperatures for earlier spectral types, but the effect of the Balmer jump makes its 
dependence with temperature non-monotonic, especially for high-gravity stars. 

A number of sensitivity curves for the U BV system have been published over the years 
and a review can be found in Bessell (1990). That author found a relatively good agreement 
between different sources for the definitions of B and V but important differences in the 
definition of U. This can be seen by comparing the sensitivity functions proposed by Bessell 
(1990) with those of Buser & Kurucz (1978). The latter authors used the definitions for 
B and V by Azusienis & Straizys (1969) and derived a new sensitivity function for U. As 
shown in Fig. 8, both sensitivity curves have similar red edges but the blue edge extends to 
considerably shorter wavelengths for the Bessell (1990) one^. 



''AH sensitivity functions in Fig. 8 and elsewhere in this article are expressed in photon-counting form 
(see Appendix A). 
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Another important feature of both the Buser Sz Kurucz (1978) and Bessell (1990) studies 
is that two definitions for the B sensitivity function are given depending on whether they are 
intended for computing U — B or B — V. Hence, each of those papers gives four different 
sensitivity functions: Buser & Kurucz (1978) names them U3, B2, B3, and V, with U — B = 
U3 - B2 and 5 - \/ = B3 - V, while Bessell (1990) names them UX, BX, B, and V, with 
U — B = UX — BX and B — V = B — V. This practice can be traced back to Azusienis & 
Straizys (1969) as an attempt to correct for atmospheric extinction but it is quite obvious 
that it is unphysical: B is B, independently of whether it is used to calculate U — B or 
B-V. 

In this section 1 test the sensitivity curves for the Johnson UBV system. In order to 
build the photometric sample I used the GCPD and followed a selection procedure similar 
to the one in the previous section for the Stromgren system. The only relevant differences 
are: 

1. I compiled one magnitude (V) and two colors {B — V and U — B). For the four white 
dwarfs in Bohlin (2000) I used the values of V in that reference and for Vega I used 
the value in Bohlin & Gilliland (2004b). The rest of the V magnitudes and all of the 
colors are from the GCPD. 

2. As opposed to Stromgren, the samples for V, B — V and U — B were built indepen- 
dently. However, the rules regarding the requirement of three data points for each star 
and the existence of 0.035 mag cutoffs in the random uncertainties were maintained. 

3. The initial values for £b-v, and Eu^b were set to zero. 

After the selection process, 108, 111, and 101 stars were present in the V, B — V and 
U — B joint NGSL + Bohlin samples, respectively. 96 stars were present in the three samples 
and the 101 stars in the U — B sample were also included in the B — V sample. The numbers 
for V, B, and B — V including only NGSL stars are 98, 104, and 96, respectively. The same 
numbers including only Bohlin stars are 10, 7, and 5, respectively. 

5.2. B-V 

The results lor B — V can be seen in Fig. 9. The plot uses the B and V sensitivity 
curves of Bessell (1990); I tried using the Buser & Kurucz (1978) and found that the results 
were very similar. £b-v was iteratively determined to be 0.012 mag. The data does not show 
any significant trend with color, with a slope of the weighted hnear fit of —0.003 ± 0.003, 



-14- 



just 1 sigma away from zero. Also, no significant difference is observed between tfie NGSL 
and Bohlin samples. Therefore, I conclude that the Bessell (1990) sensitivity curves for B 
and V provide an accurate description for those filters. 

The zero point was calculated using inverse variance weighting to obtain ZPb-v — 
0.010 ± 0.001 mag. For the systematic uncertainties I obtain 0.002 mag from the slope of 

the weighted linear fit and 0.003 mag from the temperature uncertainty in the white dwarf 
calibration. Therefore, I finally obtain ZPb-v = O.OlOiO.OOl (random) ±0.004 (systematic) 
magnitudes. It is interesting to point out that the photometric value for Vega itself compiled 
from the GCPD is B — V = 0.000 ± 0.012 (with the uncertainty coming from Eb-v), which 
is less than 1 sigma away from ZPb-v- This is an additional confirmation of the consistency 
of the results. 



5.3. U -B 

I performed for U — B an analysis similar to the one ioT B — V. Following the Stromgrcn 
mi and Ci cases, 1 chose the equivalent oi h — y (i.e. B — V) as the comparison parameter 
since it provides a quasi-monotonic function of the spectral slope in the optical range. Given 
the differences between the Buser & Kurucz (1978) and Bessell (1990) U sensitivity curves, 
the two cases were analyzed independently and the corresponding plots are shown in Fig. 10. 

The Buser & Kurucz (1978) plot (left panel of Fig. 10) shows a clear color term, with the 
bluer stars below the mean and the redder ones above it. The measured slope is 0.016±0.004, 
which is indeed 4 sigmas away from zero^°, leading to a variation over the full range of ±0.016 
magnitudes, which is too much to be treated simply as a systematic error. Therefore, I 
conclude that the Buser & Kurucz (1978) sensitivity curves do not provide an adequate 
representation of the hterature U — B data. 

The Bessell (1990) plot (right panel of Fig. 10) also shows a color term, though weaker 
and of the opposite sign as the Buser & Kurucz (1978), with the bluer stars above the mean 
and the redder ones below it. The measured slope is — 0.008 ± 0.004, which is 2 sigmas away 
from zero and yields a variation over the full range of ±0.008 magnitudes. Such a variation is 
not very large but two additional issues exist. First, as previously mentioned, the calculation 
of the U — B color using the Bessell (1990) definitions implies the unphysical use of a different 
sensitivity curve for B than the one used for the same filter in the computation of the B — V 



^°As it was the case for our analysis of Stromgren ci the exact uncertainty in the slope depends on the 
value of su-B, which is determined later. 
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color. Second, the existence of a significant slope in the weighted linear fit to the data in 
the right panel of Fig. 10 is a manifestation of a more complex effect. As we go from left to 
right in the plot, we see first that the OB stars tend to be above the mean, then the A stars 
tend to be below it, the F and G again above it, and finally the K and M stars below it. 
This indicates that the considered U sensitivity curve is not correctly assigning weights to 
the flux to the right and left of the Balmer jump. A consequence is that the derived ZP[/_b, 
0.045, is very different from the photometric U — B color for Vega, the reference star (AO V), 
which is —0.004 magnitudes. Those arguments lead me to conclude that the Bessell (1990) 
sensitivity curves do not represent the literature U — B data correctly, either. 

The inadequacy of the existing U sensitivity curves, which 1 previously discussed in 
Maiz Apellaniz (2005b), led me to derive a new one using the same procedure I followed 
for Stromgren u. In this case I applied minimization to the {U — S)phot — {U — B)spec 
data and used 12 pivot wavelengths at 100 A intervals, setting the two extremes (where 
the sensitivity goes to zero) at 3 050 A and 4 150 A, respectively. Also, I eliminated the 
unphysical definition of an alternative B filter and used the same one as for the B — V zero 
point calculation. As it was the case for Stromgren u, several iterations with different values 
of Su-B were required until a reduced of 1-0 was reached for Su-b = 0.011. 

The new U sensitivity curve is shown in Fig. 8. As it happened for the Stromgren u 
case, the red side of the curve is quite similar to the previous curves but the blue side is 
quite different. The function shows a steep rise around 3 100 A followed by a plateau and a 
new rise around 3 600 A. When the relative weights of the regions to the left and right of the 
Balmer jump are calculated, 1 find that the weight of the region to the left is smaller than for 
Bessell (1990) but larger than for Buser & Kurucz (1978). This is an expected effect, given 
the opposite slope of the weighted hnear fits of the two plots in Fig. 10: a sensitivity curve 
with near-zero slope should lie at an intermediate point between the two previous ones. It 
should be pointed out that during the fitting procedure I found out that the exact shape of 
the sensitivity curve to the left of 3 600 A is poorly constrained by the data, since alternative 
solutions (e.g. one with a small secondary peak maximum 3 300 A) had similar values of x^- 
This is also an expected result, since the 3100-3 500 A region contains little information in 
the form of strong absorption lines or edges. However, all of those alternative solutions had 
nearly constant integrated weights for the regions to the left and the right of the Balmer 
jump and, therefore, produced little variation in the derived synthetic U — B. Furthermore, 
it is important to remember that, as I mentioned in the introduction, it is likely that no 
single sensitivity curve can fit all the existing data: the purpose of this calculation is to find 
the sensitivity curve that minimizes systematic errors. 

The new version of the {B — V)phot vs. {U — -B)phot — {U — i?)spec is shown in Fig. 11. As 
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expected, the color term is absent, with the slope of the weighted hnear fit being 0.001±0.004. 
The value for ZPu_b, obtained from the sensitivity-curve minimization procedure, is 
0.020 ± 0.006 magnitudes. The systematic uncertainty from the fitted slope is very small 
(< 0.001 mag) and the one from the uncertainty in the white dwarf temperature scale is 
0.007 magnitudes. ZP[/_b is now considerably closer to the measured photometric U — B 
color of Vega but, for £u-b = 0.011, is still 2 sigmas away. Also, it should be noted that 
the five stars in the Bohlin sample appear to be shifted downward by ~ 0.02 magnitudes 
with respect to the mean locus of the NGSL sample in Fig. 11. It is unhkely that the source 
of the relative displacement is due to a problem with an incorrect relative calibration of 
spectrophotometry of the NGSL sample with respect to the Bohlin sample (as it was the 
case for Tycho-2 in Mai'z Apellaniz 2005a) since such a displacement should also show up in 
Fig. 7, given the similar wavelengths sampled by U and u. The effect could simply be a case 
of small number statistics (the probability that all of five points in a Gaussian distribution lie 
at the same side of the mean is 1/16) but, to be conservative, 1 suggest that the systematic 
uncertainty be doubled to take into account the possibility that the difference between the 
results for the NGSL and the Bohhn samples is real. Therefore, the proposed value for 
ZPu-B is 0.020 ± 0.006 (random) ±0.014 (systematic) magnitudes. 



5.4. V 

Up until now I have dealt with the calibration of colors and indices but not of magnitudes 
themselves. In the first case the accuracy of the calibration depends ultimately on the 
parameters (fundamentally temperature) of the spectrophotometric standards and on the 
stability and knowledge of the observational photometric conditions. The same is true for 
magnitudes but there one also runs into the spectrophotometric absolute flux calibration. 
That issue has been tackled by Bohlin (2000) and Bohlin & Gilliland (2004b), who combined 
the 5 556 A absolute fiux of Vega of Megessier (1995), Landolt V photometry for the three 
white dwarfs used as primary calibrators, and STIS spectrophotometry to derive a V value 
for Vega of 0.026 ± 0.008 magnitudes. 

Trying to use the Bohlin sample to derive a value for ZPy using the same procedure 
as for e.g. B — V would be falling into a logical inconsistency: the V photometry of the 
white dwarf calibrators and their HST spectrophotometry are not independent, since the 
absolute calibration of the second has been derived using the first as input. Therefore, the 
fact that the Iphot — Vspec values for the Bohlin sample in Fig. 12 show no trend with Vphot or 
{B — y)phot is simply a verification that the spectrophotometric calibration is self-consistent. 
A similar argument takes place with the value of ZPy. In the Bohlin sample there are five 
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stars where the V magnitudes are taken from the two references above plus another four 
stars with accurate photometry in the GCPD (however, only six stars have also accurate 
B — V photometric colors in our data). The weighted mean of V^phot — Vspec yields 0.024, 
which is very close to the V magnitude of Vcga^^. And what about using the NGSL sample 
to calculate ZPy? Unfortunately, there one runs into the same problem as for ZPy^ or 
ZPsrj. in Paper I: the NGSL data cannot be used to directly calculate the zero point of a 
magnitude because they were obtained with the 52X0.2 slit, so significant hght losses due to 
miscenterings are to be expected. Indeed, that effect is clearly seen in Fig. 12 in a similar 
manner to what happens in Fig. 4 in Paper I. Therefore, the most appropriate choice for 
ZPv is the V magnitude for Vega, 0.026 ± 0.008, derived by Bohhn & GiUiland (2004b). 



5.5. A new look at the absolute calibration of Tycho-2 photometry 

In Paper I the Vr filter was calibrated through the photometry for the Bohlin sample 
and the ZPbt-Vt derived from the combination of the NGSL and Bohhn samples. Since the 
last term is now known to have been in error, a new value for ZPy^ needs to be calculated. 
If I use the four stars in the Bohlin sample with accurate Vt photometry and calculate the 
mean of Vrphot ~ ^spec using inverse variance weights I obtain 0.039 ± 0.008 magnitudes 
(see right panel of Fig. 4 in Paper 1). Note that the uncertainty originates in the Tycho- 
2 photometry, not in the the Vega V magnitude. It is interesting to notice that if the 
absolute HST calibration were to be off by a multiplicative constant (let's say 1%), it would 
not affect that measurement of ZPy^ because the zero point is measured by comparing two 
spectrophotometric fluxes on the same scale (e.g. the multiplicative constant would affect 
both the numerator and the denominator in Eq. Al)^^. A different issue, of course, would 
be the absolute zero point as defined by Eq. B2. 

The availability of the Johnson V photometry allows me to perform an independent 
measurement of ZPy^. Vr and V measure the flux in a similar wavelength range (Bessell 2000) 
but are also different enough that Vr — V varies by 0.26 magnitudes as we move along the 
unreddened main sequence. In any case, since the Vt and V photometric measurements are 
independent, I can compile photometric Vr — V values and compare them with the computed 
spectrophotometric ones to derive ZPy^.y- The results are shown in Fig. 13. There is a 
color term present but the slope is quite small, 0.008 ±0.004, and is only 2 sigmas away from 



^^The reason why it is not identical is the inclusion of the five additional stars with GCPD photometry. 

^^This is the same reason why the 0.7% uncertainty in the absolute flux of Vega (Megessier 1995) is not 
included in the error balance for ZPy. 
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zero so, as I have done for other colors, the effect can be simply included as a systematic 
uncertainty. Using inverse variance weights I obtain ZPvt^-v = 0.000 ± 0.002 magnitudes. 
For the systematic uncertainty there are 0.005 magnitudes due to the color term and 0.001 
magnitudes due to the uncertainty in the white dwarf temperature scale, leading to a final 
value of ZPvt-v = 0.000 ± 0.002 (random) ±0.005 (systematic) magnitudes. ZPy^ can now 
be obtained by adding the ZPy value from the previous subsection. Merging random and 
systematic uncertainties for simplicity, I get ZPy^ = 0.026 ± 0.010. 

The measurements of ZPy^ in the two previous paragraph are independent and compat- 
ible, as they arc just one sigma away ([0.039 - 0.026]/ VO.OOS^ + 0.010^ = 1.0). Furthermore, 
since there are no systematic uncertainties in common in the two derivations, they can be 
combined using inverse variance weights to obtain a best value for ZPvt of 0.034 ± 0.006. 

The value for ZP^^ derived in Paper 1 from the Bohlin sample is 0.078 ± 0.009. If I 
combine this with the new value for ZPy^ (derived from the NGSL Vr — V colors and from 
Bohlin & GiUiland 2004b) I obtain an alternative measurement for ZPbt,-Vt of 0.044 ±0.011 
magnitudes, which is consistent with the result of 0.033±0.001 (random) ±0.005 (systematic) 
previously derived, thus providing an additional check on the validity of the zero points. 

6. Summciry and applications 

I have analyzed the sensitivity curves for the Tycho-2, Stromgren, and Johnson photo- 
metric systems and have determined that for 7 of the 9 curves there are values in the hterature 
consistent with the available photometry and HST spectrophotometry. The recommended 
photon-counting sensitivity curves are given in Tables 1,2, and 3. The results for Bt and Vt 
are taken from Bessell (2000), for Stromgren vby from Matsushima (1969), and for Johnson 
BV from Bessell (1990) and have been converted in all cases into photon counting curves. 
The results for Stromgren u and Johnson U have been derived in this paper after finding 
that the published curves do not provide consistent results with both the photometry and 
the spectrophotometry. For those two cases, the differences between the literature and the 
new curves reside on the blue edge of a filter centered around 3400-3 700 A, which indicates 
that the likely culprit is an incorrect characterization of atmospheric extinction. 

I show in Tables 4 and 5 the recommended zero points for the colors, indices, and 
magnitudes calculated in this article. All zero points use Vega as the reference spectrum. In 
Table 5 I also include the zero points for y computed by Holberg & Bergeron (2006)^^ and 



1 could not compute ZPj,-y from the data in this paper because the GCPD provides Johnson V magni- 
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for the 2MASS system computed by Cohen et al. (2003). 

Recently, Holberg & Bergeron (2006) have combined optical/IR photometry and SED 
models of very low extinction DA white dwarfs with the Bohlin (2000) HST flux calibration 
to derive results similar to the ones in this paper. Their color/index zero points for the 
Stromgren and Johnson systems are independent from the ones here because they are based 
on a largely different sample and because they use model SEDs instead of observed ones. 
Their values for ZPi,_y, ZP^^, and ZPb-v are 0.004, 0.154, and 0.007 magnitudes, respec- 
tively, all of which agree within 1 sigma with the ones in this paper. Since those authors do 
not attempt to fit Stromgren u or Johnson U sensitivity curves, it is not possible to compare 
their zero points for Ci and U — B directly with mine. However, the fact that their Stromgren 
data show the highest scatter for the u band and that their Johnson filter with the largest 
dispersion is U points in the same direction as my conclusion that new sensitivity curves are 
required in those two cases. Therefore, the results in the two articles are consistent. 

The most direct application of this work is the use of the zero points in Tables 4 and 5 
and of the two new sensitivity curves for synthetic photometry. The latest version of CHORI- 
ZOS (Maiz Apellaniz 2004) at the time of this publication (2.0) incorporates them. Another 
application is to provide uncertainties for the Stromgren and Johnson colors and indices 
in the GCPD. To do the latter, I start with the photometry originally compiled in sec- 
tions 4 and 5. For each data point for each star in the NGSL sample I can now compute the 
difference between the photometric and spectrophotometric values for each color and index 
and obtain the corresponding distributions. The total number of data points for each of the 
five colors/indices turns out to be between 662 and 871, much larger than the numbers used 
to test the sensitivity curves because now all stars in the sample are used and no average 
is done for each star; in this way one can obtain an estimate for the typical uncertainty 
associated with a hterature value. The distributions are shown in Figs. 14 and 15 and can 
be well approximated by a Gaussian plus tails. The tails are defined as being more than 3.5 
sigmas away from the mean, with the standard deviation calculated iteratively. The value 
of 3.5 sigmas is chosen because 1/[1 — errorf(3.5/\/2)] = 2149, implying that a sample of 
roughly three times the ones we have here are required for at least one data point to be at 
larger distances from the mean. For each color/index I find that 2-6% of the data points are 
in the tails and I take those objects to be either misidentifications or cases with incorrect 
data reductions. Hence, they are not considered for the subsequent analysis. 

The mean for the distributions of (6 - t/)phot - {b - |/)spec, "^iphot - "^ispec, and Ciphot - 
Cispec are all within 0.001 magnitudes of the corresponding zero points. The dispersions. 



tudes instead of Stromgren y magnitudes and a small color term y — V exists, just as for Vt — V. 



-20- 



which can be taken to be the estimates for the uncertainty associated with a single data point 
in the GCPD, are 0.013, 0.014, and 0.018 magnitudes, respectively. These results provide 
an additional check to the validity of the sensitivity curves and zero points calculated in this 
article^^ and show that the typical precision for a single value of a Stromgren color or index 
in the literature is 1 — 2%. In Tables 6, 7, and 8 the dispersions are given as a function of V 
magnitude and N, the number of measurements used for an individual data point (Hauck & 
MermiUiod 1998). The variation with V and N is in the expected direction (brighter stars or 
those with more measurements have lower dispersions) but the effect is very small, so using 
a single value for the photometric uncertainty should not be a bad approximation. 

The mean for the distributions of {B — V)phot~(-B — V")spec and {U — -B)phot — (f^ — -B)spec 
are 0.0016 and 0.0022 magnitudes larger than the respective zero points, shghtly larger than 
for the Stromgren cases but still within the values expected from the uncertainties derived 
for ZPb-v and ZP[/_b. The associated dispersions are 0.020 and 0.028 magnitudes'^, also 
slightly larger than in the Stromgren cases. The results in Tables 9 and 10 show a weak 
dependency of the dispersions with V and a negligible dependence with N, so a single value 
for the photometric uncertainty should also be appropriate for Johnson photometry. 1 con- 
clude that the typical precision for B — V and U — B values in the literature is in the 2 — 3% 
range, slightly larger than for Stromgren. 

I would like to thank Ralph Bohlin and an anonymous referee for useful comments on 
this paper, Jay Holberg for granting me access to his results prior to publication and for 
useful suggestions, Steve Willner for bringing into my attention an issue with the 2MASS 
zero points, and Charles Profhtt and the rest of the STIS team at STScI for their help with 
the spectrophotometric calibration. 
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When computing synthetic magnitudes from spectrophotometry, one has to be careful 
with an aspect that has generated some confusion in the past. For a photon- counting detector, 
such as a CCD, given the total-system dimensionless sensitivity function -P(A), the SED of 
the object /a(A), the SED of the reference spectrum /A,ref(''^)) the zero point ZPp for 
filter P, the corresponding magnitude is given by: 



The A inside the integrals in Eq. Al is present because /a (A) is in units of energy per 
time per unit area per wavelength and one needs to multiply by X/{hc) in order to convert 
from energy to photons (the he disappears because it is a constant that appears on both 
integrals). For an energy- counting (or energy-integrating) detector, the A is not present and 
the corresponding magnitude is given by: 



A. Photon- and energy-counting detectors 




(Al) 




(A2) 



This preprint was prepared with the A AS IM^jX macros v5.2. 
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where P'{X) and ZPp, are the sensitivity function and the zero point, respectively. 

It is easy to verify that one can use either Eq. Al or Eq. A2 independently of the 
physical type of detector involved if the sensitivity function is changed accordingly. For 
example, suppose I want to use Eq. Al with an energy-counting detector. If one multiplies 
and divides by A the quantities inside the two integrals in Eq. A2 and defines -P( A) = P'{X) / A 
and ZPp = ZPp, to be the photon-equivalent sensitivity function and zero point, respectively, 
then Eq. Al is recovered (this is the reason why previous works, such as Bessell et al. 1998, 
can be correct even though they use energy-integrating sensitivities). 

The problem that has arisen sometimes is that Eq. Al has been used with an energy- 
counting sensitivity function without first dividing it by A to convert it to a photon-equivalent 
sensitivity function. The likely origin of the confusion has been the old practice of defining 
sensitivity functions for photomultipliers as energy-counting combined with the present ex- 
tensive use of "black-box" synthetic photometry codes such as synphot designed for photon- 
counting detectors. 



B. Absolute and relative zero points 

ZPp, the zero point used in the previous section is a relative zero point, since a reference 
spectrum is used to define it. In a number of magnitude systems (typically, those defined 
with synthetic photometry in mind), such as VEGAMAG, STMAG, or ABMAG (STScI 
1998), ZPp = 0.0 by definition and the calibration information is contained in the reference 
spectrum alone. For other magnitude systems derived from standard stars, ZPp is not 
exactly zero and has to be measured (that is one of the main purposes of this paper), 
although it is possible that if the right reference spectrum is used (e.g. Vega) ZPp can be 
indeed approximately zero. However, one can also avoid altogether the reference spectrum 
in the definition of a magnitude system by using (for a photon-counting detector): 

mp = -2.51ogio y PiX)fx{X)XdX^ + AZPp, (Bl) 
where AZPp is the absolute zero point for filter P, which is related to ZPp by: 

AZPp = 2.5 logio Q P(A)/vef (A) A rfA^ + ZPp. (B2) 



Note that the quantities inside the parenthesis in Eqs. Bl and B2 are not dimensionless 
(they have dimensions of energy per time per length), so the units used must be specified 



when values for AZPp are given. 
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Fig. 1. — Comparison between photometric and spectrophotometric i?T — Vr colors as a 
function of photometric i?T — Vt for the NGSL + Bohhn samples. The error bars represent 
the photometric uncertainties and the horizontal blue line marks the proposed ZPs^_y^. 
The green line shows the result of a weighted linear fit to the data. 
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Fig. 2. — Histograms for the comparison between photometric and spectrophotometric 
-Bt — colors for the NGSL sample, (left) Regular histogram, (right) Histogram for the 
data shifted by the proposed ZP^^.y^ and normalized by the individual uncertainties. A 
Gaussian distribution with /i — and o" = 1 is overplotted for comparison. 
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Fig. 3. — Comparison between photometric and spectrophotometric b — y colors as a func- 
tion of photometric b — y for the NGSL + Bohhn samples. The error bars represent the 
photometric uncertainties and the horizontal blue line marks the proposed ZP b-y. The green 
line shows the result of a weighted linear fit to the data. 
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Fig. 4. — Comparison between photometric and spectrophotometric mi indices as a function 
of photometric h — y for the NGSL + Bohhn samples. The error bars represent the photo- 
metric uncertainties and the horizontal blue line marks the proposed ZP^i- The green line 
shows the result of a weighted linear fit to the data. 
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Fig. 5. — Comparison between photometric and spectrophotometric ci indices as a function 
of photometric b — y for the NGSL + Bohhn samples using the Matsushima (1969) sensitivity 
curves. The error bars represent the photometric uncertainties and the horizontal blue line 
marks the weighted mean for the vertical coordinate. The green line shows the result of a 
weighted linear fit to the data. 
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Fig. 6. — Matsushima (1969) photon-counting normalized sensitivity curves for the 
Stromgren u, v, and b filters and the new curve for u derived in this article. Three nor- 
malized SEDs from the two samples used to calibrate the Stromgren photometry are also 
shown as fx- 
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Fig. 7. — Comparison between photometric and spectrophotometric ci indices as a function 
of photometric b — y for the NGSL + Bohhn samples using the sensitivity curves proposed 
in this article. The error bars represent the photometric uncertainties and the horizontal 
blue line marks the proposed ZPc^. The green line shows the result of a weighted linear fit 
to the data. 
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Fig. 8.— Bessell (1990) U (UX) and B (B) and Buser & Kurucz (1978) U (U3) photon- 
counting normalized sensitivity curves and the new curve for U derived in this article. Three 
normalized SEDs from the two samples used to calibrate the Johnson photometry are also 
shown as fx. 
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Fig. 9. — Comparison between photometric and spectrophotometric B — V colors as a func- 
tion of photometric B — V for the NGSL and Bohhn B — V samples. The error bars represent 
the photometric uncertainties and the horizontal blue line marks the proposed ZPb-v- The 
green line shows the result of a weighted linear fit to the data. 
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Fig. 10. — Comparison between photometric and spectrophotometric U — B colors as a 
function of photometric B — V for the NGSL and Bohhn U — B samples. The left panel 
uses the sensitivity curve definitions of Buser & Kurucz (1978) and the right panel those of 
Bessell (1990). The error bars represent the photometric uncertainties and the horizontal 
blue line marks the weighted mean for the vertical coordinate in each panel. The green line 
shows the result of a weighted linear fit to the data in each panel. 
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Fig. 11. — Comparison between photometric and spectrophotometric U — B colors as a func- 
tion of photometric B ~V for the NGSL and Bohhn U — B samples using the sensitivity 
curve definitions proposed in this article. The error bars represent the photometric uncer- 
tainties and the horizontal blue line marks the proposed ZPu_b- The green line shows the 
result of a weighted linear fit to the data. 
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Fig. 12. — Comparison between photometric and spectrophotometric V magnitudes as a 
function of the photometric values for V (left) and B — V (right) for the NGSL -|- Bohlin 
samples. The error bars represent the photometric uncertainties and the horizontal line 
marks the proposed ZPy. The asymmetric scatter of values below the line is due to light 
loss at the slit due to poor centering. 
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Fig. 13. — Comparison between photometric and spectrophotometric Vr — V magnitudes as 
a function of the photometric values i?T — Vr (right) for the NGSL sample. The error bars 
represent the photometric uncertainties and the horizontal line marks the proposed ZPy^_y. 
The green line shows the result of a weighted linear fit to the data. 
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Fig. 14. — Comparison between the photometric and spectrophotometric Stromgren values 
for the individual data points in the GCPD. The left, center, and right panels show b — y mi, 
and Ci, respectively. Gaussians with zero mean and a of 0.013, 0.014, and 0.018 magnitudes, 
respectively, are plotted for comparison. 
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Fig. 15. — Comparison between the photometric and spectrophotometric Johnson values for 
the individual data points in the GCPD. The left and right panels show B — V and U — B, 
respectively. Gaussians with zero mean and a of 0.020 and 0.028 magnitudes, respectively, 
are plotted for comparison. 
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Table 1. Recommended photon-counting sensitivity curves for the Tycho B^Vr system 



A (A) 


P(X) 


A (A) 




3500 


0.000 


4550 


0.000 


3550 


0.014 


4600 


0.022 


3600 


0.058 


4650 


0.115 


3650 


0.123 


4700 


0.302 


3700 


0.206 


4750 


0.532 


3750 


0.305 


4800 


0.740 


3800 


0.416 


4850 


0.873 


3850 


0.530 


4900 


0.944 


3900 


0.636 


4950 


0.977 


3950 


0.724 


5000 


0.994 


4000 


0.787 


5050 


1.000 


4050 


0.830 


5100 


0.995 


4100 


0.861 


5150 


0.979 


4150 


0.889 


5200 


0.953 


4200 


0.920 


5250 


0.920 


4250 


0.953 


5300 


0.882 


4300 


0.982 


5350 


0.840 


4350 


1.000 


5400 


0.797 


4400 


0.976 


5450 


0.752 


4450 


0.861 


5500 


0.707 


4500 


0.685 


5550 


0.661 


4550 


0.489 


5600 


0.614 


4600 


0.317 


5650 


0.567 


4650 


0.202 


5700 


0.520 


4700 


0.136 


5750 


0.473 


4750 


0.101 


5800 


0.426 


4800 


0.080 


5850 


0.381 


4850 


0.059 


5900 


0.336 


4900 


0.036 


5950 


0.294 


4950 


0.016 


6000 


0.255 
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Table 1 — Continued 



Vv 



A (A) 


P(A) 


A (A) 


P(A) 


5000 


0.003 


6050 


0.219 


5050 


0.000 


6100 


0.187 






6150 


0.160 






6200 


0.136 






6250 


0.114 






6300 


0.097 






6350 


0.082 






6400 


0.069 






6450 


0.058 






6500 


0.047 






6550 


0.038 






6600 


0.028 
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Table 2. Recommended photon-counting sensitivity curves for the Stromgren uvby 

standard system 



uvby 



A (A) 


P(A) 


A (A) 


P(A) 


A (A) 


P(A) 


A (A) 


P(A) 


3150 


0.000 


3750 


0.000 


4350 


0.000 


5150 


0.000 


3175 


0.000 


3775 


0.003 


4375 


0.009 


5175 


0.013 


3200 


0.007 


3800 


0.008 


4400 


0.021 


5200 


0.033 


3225 


0.076 


3825 


0.018 


4425 


0.035 


5225 


0.054 


3250 


0.220 


3850 


0.031 


4450 


0.051 


5250 


0.079 


3275 


0.424 


3875 


0.048 


4475 


0.079 


5275 


0.136 


3300 


0.607 


3900 


0.064 


4500 


0.109 


5300 


0.199 


3325 


0.758 


3925 


0.100 


4525 


0.173 


5325 


0.296 


3350 


0.879 


3950 


0.162 


4550 


0.267 


5350 


0.452 


3375 


0.961 


3975 


0.268 


4575 


0.428 


5375 


0.637 


3400 


0.999 


4000 


0.410 


4600 


0.644 


5400 


0.793 


3425 


1.000 


4025 


0.609 


4625 


0.861 


5425 


0.879 


3450 


0.990 


4050 


0.813 


4650 


0.979 


5450 


0.926 


3475 


0.972 


4075 


0.960 


4675 


1.000 


5475 


0.964 


3500 


0.945 


4100 


1.000 


4700 


0.959 


5500 


1.000 


3525 


0.898 


4125 


0.970 


4725 


0.815 


5525 


0.982 


3550 


0.830 


4150 


0.878 


4750 


0.597 


5550 


0.872 


3575 


0.742 


4175 


0.748 


4775 


0.375 


5575 


0.672 


3600 


0.634 


4200 


0.568 


4800 


0.236 


5600 


0.479 


3625 


0.503 


4225 


0.367 


4825 


0.147 


5625 


0.302 


3650 


0.357 


4250 


0.224 


4850 


0.084 


5650 


0.193 


3675 


0.239 


4275 


0.135 


4875 


0.053 


5675 


0.127 


3700 


0.157 


4300 


0.080 


4900 


0.032 


5700 


0.080 


3725 


0.107 


4325 


0.053 


4925 


0.025 


5725 


0.048 


3750 


0.068 


4350 


0.039 


4950 


0.019 


5750 


0.026 


3775 


0.039 


4375 


0.027 


4975 


0.015 


5775 


0.018 


3800 


0.019 


4400 


0.014 


5000 


0.008 


5800 


0.015 


3825 


0.000 


4425 


0.007 


5025 


0.004 


5825 


0.008 






4450 


0.000 


5050 


0.000 


5850 


0.000 
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Table 3. Recommended photon-counting sensitivity curves for the Johnson UBV 



standard system 




u 


B 






V 


A (A) 


P(X) 


X (A) 


P(\) 


A (A) 


P(\) 


3050 


0.000 


3600 


0.000 


4700 


0.000 


3100 


0.237 


3650 


0.011 


4750 


0.004 


3150 


0.403 


3700 


0.033 


4800 


0.032 


3200 


0.489 


3750 


0.058 


4850 


0.084 


3250 


0.504 


3800 


0.144 


4900 


0.172 


3300 


0.508 


3850 


0.348 


4950 


0.310 


3350 


0.511 


3900 


0.601 


5000 


0.478 


3400 


0.513 


3950 


0.817 


5050 


0.650 


3450 


0.516 


4000 


0.958 


5100 


0.802 


3500 


0.528 


4050 


1.000 


5150 


0.913 


3550 


0.603 


4100 


0.995 


5200 


0.978 


3600 


0.741 


4150 


0.995 


5250 


1.000 


3650 


0.889 


4200 


0.994 


5300 


0.994 


3700 


0.985 


4250 


0.976 


5350 


0.977 


3750 


1.000 


4300 


0.950 


5400 


0.950 


3800 


0.965 


4350 


0.921 


5450 


0.911 


3850 


0.841 


4400 


0.888 


5500 


0.862 


3900 


0.648 


4450 


0.845 


5550 


0.806 


3950 


0.424 


4500 


0.792 


5600 


0.747 


4000 


0.231 


4550 


0.733 


5650 


0.690 


4050 


0.109 


4600 


0.673 


5700 


0.634 


4100 


0.035 


4650 


0.619 


5750 


0.579 


4150 


0.000 


4700 


0.569 


5800 


0.523 






4750 


0.519 


5850 


0.467 






4800 


0.467 


5900 


0.413 






4850 


0.414 


5950 


0.363 






4900 


0.362 


6000 


0.317 






4950 


0.315 


6050 


0.274 






5000 


0.272 


6100 


0.234 






5050 


0.232 


6150 


0.200 
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Table 3— 


—Continued 






TT 
U 


JD 




T/ 
V 


A (Aj I^[A) A [A) 




A {A) 




5100 


0.193 


6200 


0.168 


5150 


0.155 


6250 


0.140 


5200 


0.121 


6300 


0.114 


5250 


0.096 


6350 


0.089 


5300 


0.075 


6400 


0.067 


5350 


0.054 


6450 


0.050 


5400 


0.034 


6500 


0.037 


5450 


0.018 


6550 


0.027 


5500 


0.007 


6600 


0.020 


5550 


0.001 


6650 


0.016 


5600 


0.000 


6700 


0.013 






6750 


0.012 






6800 


0.010 






6850 


0.009 






6900 


0.007 






6950 


0.004 






7000 


0.000 
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Table 4. Color/index zero points and associated uncertainties/errors. 





Tycho-2 


Stromgren 


Johnson 


— Vt 


b-y mi ci 


B-V U-B 


zero point 


0.033 


0.007 0.154 1.092 


0.010 0.020 


random 


0.001 


0.001 0.001 0.002 


0.001 0.006 


systematic 


0.005 


0.003 0.003 0.004 


0.004 0.014 



Table 5. Magnitude zero points and associated uncertainties. The values for y are from 
Holberg & Bergeron (2006) and those for the 2MASS filters from Cohen et al. (2003). 





'iycho-2 


Stromgren 


Johnson 




2MASS 






Vt 


y 


V 


J 


H 




zero point 
uncertainty 


0.034 
0.006 


0.014 
0.010 


0.026 
0.008 


-0.001 
0.005 


0.019 
0.007 


-0.017 
0.005 



Table 6. Dispersions for {b — y)phot — {b — y)spec as a function of V and N for all the data 

points (excluding 3.5 sigma outliers). 



N 

<"3 >4 AU 



< 6.0 0.012 0.011 0.012 
V>6.0 0.014 0.013 0.014 
All 0.014 0.013 0.013 
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Table 7. Dispersions for mi p^ot — "^i,spec a function of V and N for all the data points 

(excluding 3.5 sigma outliers). 







N 




< 3 


> 4 


All 


V < 6.0 


0.013 


0.012 


0.012 


V > 6.0 


0.015 


0.013 


0.014 


All 


0.014 


0.013 


0.014 
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Table 8. Dispersions for ci p^ot — ci spec as a function of V and N for all the data points 

(excluding 3.5 sigma outliers). 







N 




< 3 


> 4 


All 


V < 6.0 


0.018 


0.015 


0.017 


V > 6.0 


0.019 


0.018 


0.018 


All 


0.018 


0.017 


0.018 
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Table 9. Dispersions for {B — V)phot — {B — V)spec as a function of V and N for all the 

data points (excluding 3.5 sigma outliers). 







N 




< 3 


> 4 


All 


V < 6.0 


0.017 


0.018 


0.018 


V > 6.0 


0.020 


0.024 


0.022 


All 


0.019 


0.020 


0.020 
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Table 10. Dispersions for {U — -B)^phot — {U — B)^spec as a function of V and N for all the 

data points (excluding 3.5 sigma outliers). 







N 




< 3 


> 4 


All 


V < 6.0 


0.026 


0.029 


0.028 


V > 6.0 


0.030 


0.027 


0.029 


All 


0.028 


0.028 


0.028 



